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Balian-Veneroni variational principle 



O 

(N 

D 



Cedric Simenel^' ^'B 
''Department of Nuclear Physics, Research School of Physics and Engineering, 
Australian National University, Canberra, Australian Capital Territory 0200, Australia 
^CEA, Centre de Saclay, IRFU/Service de Physigue Nucleaire, F-91191 Cif-sur-Yvette, France. 

(Dated: February 22, 2011) 

The Balian-Veneroni (BV) variational principle, which optimizes the evolution of the state ac- 
cording to the relevant observable in a given variational space, is used at the mean-field level to 
determine the particle number fluctuations in fragments of many-body systems. For fermions, 
the numerical evaluation of such fluctuations requires the use of a time-dependent Hartree-Fock 
(TDHF) code. Proton, neutron and total nucleon number fluctuations in fragments produced in 
collisions of two *''Ca are computed for a large range of angular momenta at a center of mass energy 
Ec.m. ~ 128 MeV, well above the fusion barrier. For deep-inelastic collisions, the fluctuations cal- 
culated from the BV variational principle are much larger than standard TDHF results, and closer 
to mass and charge experimental fluctuations. For the first time, correlations between proton and 
neutron numbers are determined within a quantum microscopic approach. These correlations are 
shown to be larger with exotic systems where charge equilibration occurs. 
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The quantum many-body problem is the root of many 
theoretical fields aiming at describing interacting parti- 
cles such as electrons in metals, molecules, atomic clus- 
ters, Bose- Einstein condensates, or atomic nuclei 
However, it can be solved exactly for simple cases only. 
The Balian-Veneroni (BV) variational principle J2] offers 
an elegant starting point to build approximations of the 
many-body dynamics and has been applied to different 
problems in nuclear physics [§-H|, hot Fermi gas ,9,], 0'* 
theory [l^, and Boson systems [ill, [13 • In particular, 
applications to deep-inelastic collisions (DIG) of atomic 
nuclei should be of interests to the upcoming exotic beam 
facilities. DIG will be used to investigate the role of the 
isospin degree of freedom in reactions and to extract the 
density dependence of the symmetry energy. Such reac- 
tions will produce nuclei in extreme, sometimes unknown, 
states (e.g., rotating nuclei with a neutron skin, or nuclei 
at or beyond the drip-lines). 

Assuming an initial density matrix Dq at to, The BV 
variational principle optimizes the expectation value of 
an observable (Q) = Tr{DQ) at a later time ti. In this 
approach, both the state D{t) and the observable Q{t) 
vary between to and ti within their respective variational 
spaces. In most practical applications, mean-field models 
are considered in a first approximation, and, eventually, 
serve as a basis for beyond- mean-field approaches [l3[ . [l3 |. 
For instance, restricting the variational space of D(t) to 
pure independent particle states, and the one of Qit) to 
one-body operators, leads to the TDHF equation [15| 
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where p is the one-body density-matrix and h[p] is the 
Hartree-Fock (HF) single-particle Hamiltonian. Accord- 
ing to this variational approach, TDHF is, then, the best 
mean-field theory to describe expectation values of one- 



body observables. However, it should not be used, in 
principle, to determine their fluctuations and correlations 



crxY = \/{XY) - {X){Y), {X and Y are one-body op- 
erators, and fluctuations correspond to the case X = Y) 
because the XY term is outside the variational space of 
the observable. Indeed, the TDHF expression for (Jxy, 



a^xvih) =Ti-{Yp{h)X [I - p{h)]}, 



(2) 



where / is the identity matrix, has been tested on frag- 
ment mass and charge fluctuations in DIG 16|. In this 



case, X — Y counts the nucleons or protons of one frag- 
ment in the exit channel. X and Y are the matrices asso- 
ciated to X and Y, respectively, in single particle space. 
It was shown that TDHF strongly underestimates exper- 
imental fluctuations Il6i. This is an intrinsic limitation 



to TDHF [17[ which can be understood by the fact that, 
in such violent collisions, transfer of many particles may 
occur leading to very different mean-fields than the quasi- 
elastic one, while TDHF assumes that all mean-fields are 
the same. As a result, transfer of many particles are arti- 
ficially hindered. The knowledge of such fluctuations is, 
however, crucial to all quantum systems. Thus, their the- 
oretical prediction is an important challenge for quantum 
many-body models. 

To optimize fluctuations of one-body operators, the 
variational space for the observable has to be increased 
to Q e {e^"*"}, where 7 G IR and d and a'^ are parti- 
cle annihilators and creators, respectively 18, 1^. This 
leads to a prescription 



2 N r Tr{[p(to) ~ px(to,e)] [p{ta) - py(^o,e)]} 
axYih) = hm — 

(3) 

differing from Eq. ([2]). The one-body density matrices 
px{t, e) obey the TDHF equation ([T]) with the boundary 
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condition 



px{ti,t) = exp(ieX)p(ii)exp(-ieX), 



(4) 



while p{t) is the TDHF solution with initial condition 
p{tQ) being the one-body density matrix of Dq. The re- 
sult in Eq. ^ takes into account possible fluctuations 
around the TDHF mean-field evolution in the small am- 
plitude limit, i.e., at the RPA level 18, [3| [see also 
Ref. [2^ for an alternative derivation of Eq. Q]. These 
fluctuations are generated by the boost in Eq. (HJ and 
propagated in the backward Heisenberg picture from ti 
to to according to the dual of the time-dependent RPA 
equation. This is why <TxY{ii) is expressed as a func- 
tion of density matrices at the initial time t^. It is easy 
to show that, if the backward trajectories px have the 
same mean-field as the forward evolution, then Eq. ([3]) 
leads to the TDHF expression in Eq. ©. If, however, 
(small) deviations occur around the original mean-field, 
then additional terms appear and lead to an increase of 
cxviti)- In the following, TDHF and BV fiuctuations 
or correlations refer to Eq. ([2]) and ([3]), respectively. 

In this work, the BV variational principle is used in 
realistic calculations of heavy-ion collisions and first de- 
tailed comparisons with experiments are performed. The 
technique is similar to the one employed in 0, 0| • The 
neutron, proton and mass fluctuations (Tattv, czz, and 
aAA, respectively, and the correlation gnz between neu- 
tron and proton numbers distributions, are computed in 
fragments resulting both from deep-inelastic and quasi- 
elastic collisions. The correlations a^z are determined 
for the first time within a quantum microscopic approach. 

The tdhfSd code is used with the SLy4(i parameteri- 
zation 2l| of the Skyrme energy-density- functional [22 1. 
The TDHF equation ([T]) is solved iteratively, with a time 
step At = 1.5 X lO"^** s, in the center of mass (cm.) 
frame. At — 10^^** s is also used to confirm the conver- 
gence of the BV fluctuations. The single-particle wave- 
functions are evolved on a Cartesian grid of 56 x 56 x 28/2 
points with a plane of symmetry (the collision plane) and 
a mesh-size Ax = 0.8 fm. The initial distance between 
collision partners is 22.4 fm. Refs. [13, 23 1 give more 
details of the TDHF calculations. 

To evaluate Eq. ([3]) , the first step is to perform a TDHF 
evolution from tg to ti. To account for the transforma- 
tion of Eq. dl]), at time ti, the occupied single particle 
wave functions are boosted according to \ipxj{ti,(-)) = 
exp(ie(7Xj0(f))|v5j(ti)), where X stands for N, Z, or A. 
If the occupied single particle wave-function ipj refers to 
a proton (resp. a neutron), q^q. = and qz^ = 1 (resp. 
qN- = 1 and qz =0), while = 1 for protons and neu- 
trons. The function 0(r) is equal to 1 for the fragment 
on which the fluctuations are calculated, and elsewhere. 
The time ti is determined, for each collision, by the time 
at which at least one fragment cm. reaches 11.2 fm from 
one edge of the box. It ensures a minimum separation 
distance of 22.4 fm for symmetric collisions. This value is 
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FIG. 1: (a) Scattering angle as a function of angular momen- 
tum L in units of % for *°Ca+'*°Ca at S^.m. = 128 MeV. 
n is the number of times the x (collision axis) or y axis has 
been crossed by the fragment. The solid line is to guide the 
eye. The dashed line corresponds to Rutherford trajectories. 
The inset gives examples of trajectories for n = 1 and 3. (b) 
TDHF (dashed line) and BV (circles) fluctuations of Z, BV 
correlations between N and Z (triangles), and intrinsic an- 
gular momentum J/ of the outgoing fragments (dot-dashed 
line). The TKEL (solid line) is given in units of 10 MeV. The 
vertical dotted lines show the range of L leading to orbiting 
or capture trajectories. 



large enough to ensure a convergence of uxy with ti as 
the fragments interact only via Coulomb repulsion [l, [sj . 

The second step is to compute a backward evolu- 
tion from ti to to of each set of single particle wave 
functions (/jjsfi (t, e). Several values of 10^** < e < 
10^^ are considered to determine the limit in Eq ([3]). 
Following Refs. 0, [1], the initial density matrix p(to) 
in Eq. ([3]) is replaced by a backward evolved den- 
sity matrix p/(to,e — 0), i.e., without the transfor- 
mation in Eq. (|4]), to minimize numerical inaccura- 
cies. Note that the latter are easily controllable and 
this procedure is not necessary with a smaller time 
step At. The trace in Eq. (|3]) is then evaluated with 
Tr{[p/(to,0) -px(to,e)] [/3/(to,0) -py(to,e)]} = t^u + 

VxY-Vix-ViY, where ryxx' = 
and the sums run over occupied states. The quadratic 
evolution of the trace with e is used as a convergence 
check [5], as well as the property ry// = At, where At is 
the total number of nucleons. 

Let us first investigate the collision of two *°Ca nu- 
clei at Ec.m. = 128 MeV. This system is considered as 
a benchmark for experimental studies of DIC 24 , |2^ . 
Theoretically, the HF ground-state of ^"Ca is expected 
to be a good approximation because of its doubly-magic 
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nature. Fig. [TJa) shows the scattering angle 9c.m. as a 
function of the angular momentum L. For L > 90, the 
trajectories are close to Rutherford scattering (dashed 
line) and correspond to quasi-elastic reactions with small 
total kinetic energy loss (TKEL) as shown by the solid 
line in Fig.[IJb). At lower L, deviations from the Ruther- 
ford formula occur because of nuclear attraction, lead- 
ing to rotation of the di-nuclear system formed by the 
two fragments in contact. For instance, L < 72 de- 
fines (arbitrarily) the orbiting region where the fragments 
have crossed their incoming trajectory, i.e., with n > 4 
crossings of the x (collision axis) or y— axis [see inset in 
Fig. mja)]. Following damped events are defined by 
a TKEL> 30 MeV, corresponding to i < 82 in Fig. [II 
We see that a wide range of scattering angles may oc- 
cur for these damped events, which is a known feature 
of Die. For L < 66, capture occurs. It corresponds to 
a fusion cross section of ^ 1140 mb with the sharp cut- 
off formula [l3| , in good agreement with a fit on fusion- 
evaporation measurements at lower energies (26| . 

Fig. Wljo) shows BV and TDHF predictions of azz 
(only charge fluctuations are shown for clarity). The 
TDHF fluctuations [Eq. ^] have been determined from 
the probability distributions of A, Z and N in the frag- 
ments at time ti BV predictions from Eq. ([3]) at 
L = 71 are not shown as no numerical convergence with 
e could be obtained. Strong variations of the BV fluc- 
tuations are observed in the orbiting region, but there is 
no clear relationship with the amount of orbiting quan- 
tified by n in Fig. [lla). The BV fluctuations are much 
more important than the TDHF predictions for damped 
events. However, at large L (quasi-elastic reactions), the 
BV and TDHF fluctuations are similar. 

The evolution of azz with 9c.m. is plotted in Fig. [5] 
for d amp ed events, and compared with experimental 
data [2J). TDHF fluctuations show no angular corre- 
lation and strongly underestimate data. BV fluctuations 
increase with 9c.m. at small angles and form a plateau at 
large angles, in qualitative agreement with data. Quan- 
titatively, the experimental plateau is underestimated. 
This might be attributed to fusion-flssion events leading 
to large fluctuations and not included in the calculations. 
Indeed, the compound-nucleus flssion cross-section, esti- 
mated to be ~ 280 mb [2J] , is not negligible compared to 
the cross-section for damped binary events of ~ 570 mb 
from the present calculations. Note that this fission cross- 
section seems reasonable as it corresponds to a fusion- 
evaporation cross-section of ^ 860 mb (obtained from 
the difference between the fusion and fusion-fission cross 
sections) , which is compatible with data [26| . Due to the 
isotropic distribution of fission fragments, fusion-flssion 
would mostly affect large angles in Fig. [2] and may ac- 
count for the difference between BV predictions and data. 
As the fragments cool down by nucleon emission, their 
fluctuations might also increase Although the num- 
ber of TDHF iterations is too small to allow a full decay 
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FIG. 2: Comparison between BV (solid line) and TDHF 
(dashed line) predictions of azz for damped events (see text) 
as a function of dc.m. with data (circles) from [24| . 



of the fragments by nucleon emission, we can estimate 
their excitation energy E* ~ TKEL/2 and angular mo- 
menta J/ ~ (1/ — Lout)/2 [see Fig.[Hb)], assuming equal 
sharing, where Lout is the angular momentum between 
the fragments in the outgoing channel. Calculations us- 
ing the code PACE4 27| with a level density parameter 
A/7.5 MeV~^ show that the decay of fragments produced 
in Die has only a small effect on uzz (see also [3l). 
However, the average of the fragment charge distribution 
after decay goes from ^ ~ 19 for L c± 80 down to 18 in 
the orbiting region while experimental data give Z ~ 17 
at large angle 2J|. This might also be a signature of 
fusion-fission events. Indeed, symmetric fission leads to 
E* ~ 38.2 MeV according to the Viola systematics [isf . 
and, then, to more emission of light particles than in DIG 
where E* ~ 31 MeV in average in the orbiting region [see 
Fig.HJb)]. Beyond mean- field correlations may also affect 
these fluctuations. 

Fragment mass distributions in "^^Ca-l-^^Ca have been 
measured at lower energies, Ec.m. = 98.5 and 115.5 MeV, 
where fusion-fission can be neglected [i^. In this ex- 
periment, the fragments are associated to almost fully 
damped collisions with a sin~^ Oc.m. dependence of their 
cross-section. In the present work, such collisions occur in 
the orbiting region for which (jaa — 9.7 in average. This 
is in fair agreement with the data which give ctaa — 11. 

In addition to fluctuations, the BV correlations a^z 
have been computed [triangles in Fig. [T](b)]. These finite 
values of aj^z are at variance with the TDHF correlations 
which are strictly zero because the single particle states 
are assumed to have pure isospin. In fact, the probability 
P{N, Z) to have a fragment with Z and N, in TDHF cal- 
culations, is the product of the pro babilities P{Z)P{N) 
to have Z and iV, independently [23|. For instance, in the 
symmetric collisions studied here, the TDHF probability 
to have the N = Z ^'^S nucleus is the same as for the neu- 
tron rich '*°S. The latter should, however, be hindered by 
the symmetry energy which induces a fast charge equili- 
bration in the fragments Such effect is included in 
the BV approach which give a^z — crzz in damped col- 
lisions, while, in quasi-elastic reactions, correlations are 
negligible compared to fluctuations [see Fig. [Ijb)]. Note 
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that the fact that BV and TDHF predictions are sim- 
ilar for large L allows the use of TDHF to investigate 
quasi-elastic reactions [2^ . The present calculations also 
give cr AfAT ~ cFzz for all L, gaa — ^■'zxJzz for L > 90 
and OAA — '^o'zz in the orbiting region. This is in good 
agreement with the relation a\j^ = a\z + ^ + ^o"^ 
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and the behavior of ctnz- 

Finally, the ^0,92^^^.^902;^- systems at beam energy 
E/A = 8.5 McV are investigated to study the role of 
isospin asymmetry in the entrance channel. The calcula- 
tions are performed at L = 192.2 with ^"Kr and 222.6 
with ^^Kr corresponding to touching spheres of radii 
l.2A^/^ at closest approach. The N/Z ratio of ^"Kr and 
^"Zr are similar (1.22 and 1.25, respectively) but differ 
from the one of ^^Kr (1.56). Then, charge equilibration 
occurs only with ^^Kr where the calculations give an av- 
erage N/Z ~ 1.4 in both outgoing fragments. The BV 
prescription gives azz — 5.3, gnn — 7.1, and ctjvz — 5.7 
with ^"Kr, and azz — 4.7, gnn — 8.4, and u^z — 8.5 
with ^^Kr. Fluctuations are only slightly affected, while 
charge equilibration strongly increases correlations be- 
tween N and Z distributions (by ^ 50%). 

The BV prescription leads to fragment mass and 
charge fluctuations in better agreement with experiment 
than TDHF for violent collisions, although TDHF should 
be sufficient for quasi-elastic reactions. The predictions 
of correlations between N and Z distributions is an at- 
tractive feature which should be compared with experi- 
mental data where both mass and charge of each frag- 
ment are measured. In particular, reactions with exotic 
beams will allow strong N/Z asymmetries, and such cor- 
relations are expected to be increased by the charge equi- 
libration process. Applications to multi-nucleon transfer 
in actinide collisions could be used to predict probabil- 
ities for super- heavy element production [sol l3l| . The 
role of pairing correlations on fluctuations should be 
investigated using recent time-dependent Hartree-Fock- 
Bogoliubov codes [32-34|. Stochastic- mean- field meth- 
ods might also be applied to investigate the role of initial 
beyond-mean- field correlations on fluctuations 2^, 35 1. 
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